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Abstract 

The Bethe-Salpeter amplitudes of the bound states formed by two scalar particles 
with unequal masses are analyzed in the massive scalar particle exchange ladder 
model. The norms of the amplitudes are calculated numerically, and it is confirmed 
that the norm vanishes for the bound state corresponding to the complex eigenvalue 
of the coupling constant. The behaviour of the Bethe-Salpeter amplitudes in the 
momentum space is also investigated. 
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§ 1. Introduction 

It is expected that the homogeneous Bethe-Salpeter (BS) equation describes the 
bound state in quantum field theory.0 The BS kernel appearing in the equation is usually 
approximated by the contribution from the ladder graph only, and this is referred to as the 
ladder model. In this model, the BS equation can be regarded as an eigenvalue problem 
for the coupling constant at a given mass of the bound state. 

For the scalar-scalar ladder model with scalar particle exchange, if the mass of the 
exchanged particle is nonzero, and if the two scalar particles have unequal masses, the 
corresponding eigenvalue equation cannot be reduced to a real form, when the invariant 
bound-state mass squared lies in the physical region, that is, between zero and the two- 
particle threshold. In this case, it becomes a nontrivial problem whether the eigenvalues 
of the coupling constant are real or complex. 

Naito and NakanishiU made an argument that all eigenvalues would be real provided 
that eigenvalues had a certain analytic property. In the same year, zur LinderMI published 
a paper suggesting the reality of eigenvalues by numerical calculation, and KaufmannS 
published a paper suggesting, on the contrary, that eigenvalues could become complex for 
some mass configurations of three scalar particles and the bound state. Ida! pointed out, 
in the next year, that the analyticity of eigenvalues assumed in Ref. [ZJ is never satisfied if 
the eigenvalue becomes non-real. Based on a rapid progress in computing technique made 
some twenty years in between, a rather extensive numerical calculation was performed to 
resolve these entangled situation. It was found that complex eigenvalues appeared in the 
first and second excited states for the bound state mass squared around the pseudothresh- 
old, and it was also argued that the appearance could not be regarded as an artifact of 
the numerical scheme employed there. 

As a continuation of Ref. which exclusively dealt with the behaviour of eigenval- 
ues, this paper concerns with the analysis of BS amplitudes (eigenvectors) of the eigenvalue 
equation. In the next section, the general formalism for the BS equation is presented. In 
§3, the norm of the BS amplitude is calculated numerically. The norm of the BS ampli- 
tude corresponding to a complex eigenvalue is found, to within the numerical accuracy, 
to vanish. In § 4, the BS amplitudes of some bound states are depicted in the momentum 
space. Discussions and a further outlook are made in the final section. 



§ 2. BS equation and BS amplitude 

The BS (Bethe-Salpeter or bound state) amplitude <fi(p,P4) with mass squared s (> 
0), formed by two scalar particles with masses 1 + A and 1 — A by exchanging a scalar 
particle with mass /x obeys, in the ladder model, the following BS equation (after the 
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Wick rotation) 



[(1 + Am -a)+pi- 2i(l + A)y/Bp A \ 

x I A fd 4 v> f2-n 

[(1-Z\) 2 (1- ^+^ + 22(1-^)^4] vr 2 7 %2 + (p _ p ,)2 > ^ ^ 

with a = s/4 and ^ Z\ < 1. We shall denote, hereafter, the model specified above 
as the [l + A<=/j,=^l — A] model. We regard Eq. (2-1) as an eigenvalue equation for 
the coupling strength A, which has mass dimension +2. In this model, the two-particle 
threshold is equal to 4, while the pseudothreshold is 4A 2 . 

By introducing the polar coordinates in the four- dimensional momentum space (p,p^), 

Pi = p cos (3, p = p sin /3(sin 9 cos ip, sin 9 sin ip, cos 9) , 

(p>0, 0£f3, 9£n, 0£p<2n) (2-2) 

we can reduce Eq. (2-1) to an infinite system of one-dimensional integral equations. Since 
the three-dimensional angular momentum / is a good quantum number, we can assume 
that the BS amplitude has the form 

00 

0(p,p 4 ) = J2 N l,i ■ (smP) l C l L t)(cosP)Y lm (0^)Mp) (2-3) 

L=l 

for suitable I (= 0, 1,2,.. .) and m (— / <m The four-dimensional spherical harmon- 
ics are (siri f3) 1 C^^cos f3)Yi m (9 , (p), with associated normalization constants N L ,i- 

By changing the variable p (the magnitude of the four-dimensional momentum) to z 
defined by 



p=Wi±f, Z = P±, (*:-l-l asp: O^oo) (2-4) 
V 1 — z p 2 + 1 

Eq. (2-1) is transformed into 

g w (z) = A]T/ dz'K LV ^z,z')g v ^z'). (L = 1,1 + 1, . . .) (2-5) 
The BS amplitude 4>{p,Pa) is expressed in terms of the Ql,i{. z ) 's as 

- ^(^^iyJ'^XS^t^ri)^^' (2 ' 6) 

The derivation of the above equations and the explicit form of the integral kernel 
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matrix functions Kll',i(z, z') for I = and / = 1 are given in Ref. pjl These matrix 



functions are all real, and if the eigenvalue equation (2-5) admit a real eigenvalue A we 
can take the eigenvectors gi.i{z) (L — I, I + 1, . . .) all real. Even in this case, however, 
since imaginary unit % appears on the right hand side of Eq. (2-6), the BS amplitude 
4>(p,Pa) takes a complex value in general. 



§ 3. Norm calculation 

In Ref. ^]|, we have analyzed numerically the eigenvalue problem (2-5) for the s- 
wave (I = 0) case. We have approximated the infinite system of integral equations by a 
finite system of algebraic equations: The summation over L' from to oo in Eq. (2-5) is 
truncated at some cutoff value L c , and the integration over z' is replaced by the iV-point 
Gauss-Legendre quadrature formula. By this procedure, we can reduce Eq. (2-5) to the 
(L c + 1) • iV-th order matrix eigenvalue problem: 

L c N 

g L {zj) = A J2 K LL',o(zj,z k )w k g L '{z k ) , 

L'=0fc=l 

(L = 0,1,---,L C , .7 = 1,2,.-., AT) (3-1) 

where the points Zi,---,Zn are the Gauss-Legendre points on the interval [—1,1], and 
Wi, ■ ■ ■ ,wn are the corresponding integration weights. We have put gh{.Zj) '■= gL,o{zj). 

We have found that for the [1.6 <= 1.0 =>• 0.4] model (that is, for the case A = 0.6 
and \i — 1.0 ), which is the same as the case zur LindenS and KaufmannEH treated, 
complex eigenvalues appear in the range 0.25 < s < 2.65, where s is the bound-state mass 
squared. Main calculations were performed for iV = 45 and L c = 7. We have, in this 
paper, recalculate the first three eigenvalues (denoted as A , Ai and A 2 ) for the case s = 
0(the left edge of the threshold), s = 0.23(near before the start of complex eigenvalues), 
s = 0.27(near after the start of complex eigenvalues), s = 1.44( the pseudothreshold), 
s = 2.62(near before the end of complex eigenvalues), s = 2.67(near after the end of 
complex eigenvalues), s = 3.90(near the nonrelativistic limit), taking N = 47 and L c = 9. 
The result is summarized in Table I. 

Table I. Eigenvalues of the [1.6^=1.0=^0.4] model at several s's. 



s 


0.00 


0.23 


0.27 


1.44 


2.62 


2.67 


3.90 


Ao 


3.448 


3.342 


3.323 


2.755 


2.113 


2.083 


1.109 


Ai 


16.433 


16.459 


16.526 - 0.127i 


14.887 - 0.734i 


13.059 - 0.159* 


12.794 


8.285 


A 2 


17.345 


16.702 


16.526 + 0.127? 


14.887 + 0.734z 


13.059 + 0.159i 


13.157 


12.455 
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Before calculating the norm of the BS amplitude, we shall check to what accuracy 
the eigenvalue equation is satisfied. Denoting Eq. (3-1) symbolically by g = XKg, where 
K is 470 x 470 matrix, we calculate the relative error || g — XKg || / || g ||. As the vector 
norm || v ||, we consider two cases; the Z 2 - n orm ( || v ||:= (J2i\ v i\ 2 ) 1 ^ 2 ) an d Zoo-norm ( 
|| v || := maxj \v$\ ). Numerical values of the relative errors are almost same for the Z 2 -norm 
and Zoo-norm. The relative errors corresponding to the ground states (Ao) in Table I are 
all less than 10~ 6 . Rather large errors appear in the first excited state (corresponding to 
Ai) at s = and in the first and second excited states (corresponding to Ai and A 2 ) at 
s = 0.23. The errors are, however, of order 10~ 4 . For other states, including those having 
complex eigenvalues, the relative errors are less than 10 -5 . 

As was explained in Ref. ^ our numerical method is based on the iteration-reduction 
scheme .0 The relative error defined above is not referred to the reduced eigenvalue prob- 
lem (for the excited state), but to the original eigenvalue problem: We must restore the 
eigenvector from the reduced eigenvector according to a certain prescription.0 Since the 
eigenvalues are obtained to within 0.1% relative error, it can be said that our scheme can 
reproduce the eigenvector to the full accuracy as we can expect. 

Having thus confirmed that the eigenf unction gh{zj) satisfies the eigenvalue equation 
(3-1) accurately, we will calculate the norm of the BS amplitudes corresponding to the 
states listed in Table I. By the name "norm of the BS amplitude", we mean a quantity 
defined by 



This is the same as considered by Naito-Nakanishi (the left hand side of Eq. (2-8) in 
Ref. |2]]) and by Ida (Eq. (5-37) in Ref. [5]]). This quantity can be shown to take a real 
value. By making use of the expression for the BS amplitude (Z = case of Eq. (2-6)), 
and using the recursion relation of the Gegenbauer polynomial C\, we can perform the 
(four-dimensional) angular part of the integration in Eq. (3-2), to the result 



/ 



P J P40*(p, -Pa) [(1 + Af(l -a)+p 2 - 2i(l + A)^p 4 
x "(1 - Z\) 2 (l - a) + p 2 + 2z(l - A)yfcpA 0(p,p 4 ) • 



(3-2) 




dz^\(2- 5 Lfi )(l - A 2 )a{\ - z 2 ) + [l +z+ (l- A) 2 (l - <r)(l - z)\ 



L=0 




L=0 



x [Re g L (z)Re g L +i(z) + Im ^(^Im g L +i(z)] 
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+ / dz £ 2(1 - A 2 )<j{-l) L+ \\ - z 2 ) L+3 

x [Re g L (z)Re g L +2{z) + Im g L {z)\m g L +2(z)] ■ (3-3) 

In the above expression, Re <?l(z) (resp. Im <7l(z)) means the real (resp. imaginary) 
part of If the corresponding eigenvalue is real, it is natural to take gi{ z ) as a real 

function. In this case Im giiz) = and it can be shown that the equality <p*(p, —P4) = 
4>{p,P4) holds, that is, the real (imaginary) part of the BS amplitude (f)(p, p^) is an even 
(odd) function of p^. If the eigenvalue A is complex, gL(z)'s acquire imaginary parts, and 
4>{p,P4) and <f)*(p, —Pa) become linearly independent, and these two BS amplitudes have 
mutually complex conjugate eigenvalues A and A*. 

The approximate value of the norm can be obtained by putting g^z) — for L 
larger than the cutoff value L c (= 9), and by performing the integration according to the 
N(= 47)-point Gauss- Legendre formula. Since the norm in Eq. (3-3) depends on the 
normalization of ^(^)'s, that is, an overall multiplication by a constant, it is desirable to 
remove this arbitrariness. We divide the norm by a manifestly positive quantity obtained 
by replacing each summand in the integrands in Eq. (3-3) by its absolute value. The 
result is 



Table II. Norms of the BS amplitude in [1.6^=1.0=^0.4] model at several s's. 



s 


0.00 


0.23 


0.27 


1.44 


2.62 


2.67 


3.90 




1.000 


0.992 


0.990 


0.921 


0.761 


0.750 


0.189 


m 


0.999 


0.216 


-0.186 • 10~ 5 


0.283 • 10~ 5 


-0.324 • 10~ 5 


0.735 • 10- 1 


0.238 


ri2 


-1.000 


-0.216 


-0.186- 10- 5 


0.283 • 10~ 5 


-0.324 • 10~ 5 


-0.729- 10- 1 


-0.280 



In the numerical analysis of an eigenvalue problem, it is a general feature that the 
accuracy of eigenvectors is lower than that of eigenvalues. Before starting our norm 
calculation, we expect, therefore, that the numerical value of the norm corresponding 
to the zero-norm state will be at best of order 10~ 3 . We can thus safely say that the 
norms of the BS amplitudes corresponding to complex eigenvalues vanish. Theoretically 
it is known that if the norm of a BS amplitude does not vanish, this BS amplitude 
has a real eigenvalue. Although the converse statement is not true, we regard that our 
result amounts almost to an analytical proof that in the [1 + A<=fi=>l — A] model 
there can appear complex eigenvalues. We have applied our calculation procedure to the 
[1.0 <= 1.0 =^ 1.0] model (equal-mass scalar-scalar model) and the [1.6 <= =>- 0.4] model 
(unequal-mass Wick-Cutkosky model). The norms of the first three eigenstates at the 
same energies as in Table I are found to be all nonzero. The minimum (in absolute value) 
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is ~ 1/30, which corresponds to the (degenerate) first and second excited state at the 
pseudothreshold in the [1.6^=0=^0.4] model. This value cab be regarded well away from 
zero. 

§ 4. BS amplitude in momentum space 

We have verified in the preceding section that the approximate s-wave eigenvector 
g^\z;s), corresponding to the k-th eigenvalue A = X k (k = 0,1,2) at the bound-state 
mass squared s, satisfies the eigenvalue equation (3-1) with a relative error less than 
0.001% except for a few cases. We can thus expect that the corresponding approximate 
BS amplitude 4>k(\p\,P4, s) for k — 0, 1, and 2, 

- ^T^rM p MfcH 57? (4 ' 1) 

reflects the behaviour of the original BS amplitude in the Wick-rotated four- dimensional 
momentum space rather accurately. 

As a typical example of the real eigenvalue case, Figs. 1 (a) and (b) shows, respec- 
tively, the real and imaginary part of the BS amplitude of the first excited state (k = 1) at 
s = 0.23. The graphs are plotted on the p4-\p\ half plane, where |p| denotes the magnitude 
of the three-dimensional space part of the four- dimensional momentum. The amplitude 
is normalized so that the largest value of the magnitude of its real part becomes 1.00. 

Fig. 1 (a) Fig. 1 (b) 



The expression (4-1) takes a form suitable for plotting in the polar coordinates: From the 
numerical calculation, we get the data g^\ z j \ s ) f° r £ = ~ 9 and j = —23 ~ 23. 
The corresponding magnitude of the four-dimensional angular momentum p — pj :— 
y/(l + Zj)/(l - Zj) (cf. Eq. (2-4)) takes the values p_ 23 = 0.02532 (z_ 23 = -0.9987), 
p_ 22 = 0.05817 (z_ 22 = -0.09933), • • • , p = 1.0000 (z = 0.0000), • • • , p 22 = 17.19 (z 22 = 
0.9933), P23 = 39.50 (z 2 % = 0.9987). In spite of this, we have plotted the graph in 
the Cartesian coordinates on the P4-|p| plane, relying on the third-order interpolation- 
extrapolation method, in order to make use of several useful graphic software tools. The 
graph is cut if p^ or \p\ exceeds 1.00, since the behaviour of the amplitude is quite smooth 
in the region omitted. The curves in the P4-|p| plane are contour lines of the surface, and 
we can read off from these curves that the real (imaginary) part of the BS amplitude is 
indeed an even (odd) function of p4, as was explained in §3. 

The real and imaginary part of the BS amplitude corresponding the the complex 
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eigenvalue Ai = 14.8876 — 0.734z at s — 1.44 (pseudothreshold) are depicted in Figs. 2 
(a) and (b), respectively. 



Fig. 2 (a) Fig. 2 (b) 



In this case, the symmetry (antisymmetry) of the real (imaginary) part of the amplitude 
under the p^ inversion is broken. It seems, however, impossible to judge from Fig. 1 and 
Fig. 2 that the BS norm of the former is ~ 0.216 while that of the latter is ~ 0.0. 

As the last example, the absolute value of the ground state at s = 3.9 is plotted in 
Fig. 3. 

Fig. 3 



The support of the amplitude is concentrated near the origin, especially in the p^ variable, 
as it should be in the nonrelativistic limit. 

§ 5. Summary and discussion 

In this paper, we have analyzed numerically, as a continuation of our previous paper ,0 
in which the behaviour of eigenvalues of the [1.6 1.0=^>0.4] model in the s-wave case was 
clarified, the properties of the BS amplitude (eigenvector) . We have verified that the BS 
amplitude satisfies the truncated eigenvalue equation (3-1) rather accurately, especially 
well for the amplitude of the ground state and of the state having a complex eigenvalue. 

Based of this observation, we have calculated numerically the BS norm of the am- 
plitude (Eq. (3-2)), and found out that the norm vanishes, to within our numerical 
accuracy, for the eigenvector corresponding to the complex eigenvalue. We can thus say 
with certainty that the [l + Z\^=yU=^l — A] model admits complex eigenvalues for some 
configuration of A, /z and the bound-state mass. We have also depicted the BS amplitude 
in the momentum space for some typical cases. 

A preliminary calculation suggests that the eigenvalues are all real for the p-wave 
case of the [1.6 <^= 1.0 =^0.4] model. If so, it will be of some interest to see whether this is 
a universal feature valid for other values of A and /i, or not. The appearance of complex 
eigenvalues will be closely connected with the phenomenon of anomalous threshold in the 
triangle Feynman graph. These points will be discussed in forthcoming papers. 
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Figure 1: The BS amplitude of the first excited state (corresponding to the real eigenvalue) 
at s — 0.23, (a) real part and (b) imaginary part. 
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Figure 2: The BS amplitude of the first excited state (corresponding to the complex 
eigenvalue) at s — 1.44, (a) real part and (b) imaginary part. 
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Figure 3: The absolute value of the BS amplitude of the ground state at s — 3.9. 
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